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ABSTRACT 

We continue investigation of the hidden plane-mirror symmetry in the distribution of excursion sets in cosmic microwave 
background (CMB) temperature anisotropy maps, previously noticed in the three-year data of the Wilkinson microwave 
anisotropy probe (WMAP), using the WMAP 5 years maps. The symmetry is shown to have higher significance, 
< 1.7, for low multipoles ^ < 5, while disappearing at larger multipoles, > 3.5 for i > 10. The study of the 
sum and difference maps of temperature inhomogeneity regions, along with simulated maps, confirms its existence. The 
properties of these mirroring symmetries are compatible with those produced by the Sachs- Wolfe effect in the presence 
of an anomalously large component of horizon-size density perturbations, independent of one of the spatial coordinates, 
and/or a slab- like spatial topology of the Universe. 
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1. Introduction 

The properties of the CMB temperature anisotropy, as well 
as its polarization, are among the basic sources of informa- 
tion on cosmological parameters ([de Bernardis et al. 2QQQ| 
ISpergel et al. 2QQ7| IKomatsu et al. 2008^ . Their tiny fea- 
tures, such as the local spikes in the multipoles power 
spectrum, deviation from the statistical isotropy, and non- 
Gaussianity signatures, may all be the result of various fun- 
damental processes having occurred in the early Universe. 
Among the reported anomalies are the alignment of the 
principal (Maxwellian) vectors of low multipoles, the north- 
south power asymmetry, the southern anomalous cold spot, 
etc (see de Oliveira- Costa et al. 2004; Copi et al. 2004,2007; 
Schwarz et al. 2004; Eriksen et al. 2004,2007; Cruz et 
al. 2005; Morales & Saez 2008). In the present paper we 
continue the study of another deviation from statistical 
isotropy: the hidden partial plane-mirror symmetry in the 
distribution of CMB temperature fluctuations excursion 
sets, i.e. of one-connected pixel sets equal to and higher 
than the given temperature threshold (lower for negative 
thresholds) previously found in the WMAP 3-years temper- 
ature maps (Gurzadyan et al. 2007b ). We use the WMAP 
5-years data (|Hinshaw et al. 2008[) , not only to confirm the 
mirroring effect found in WMAP3 maps while studying the 
inhomogeneities in the distribution of the excursion sets, 
but also to reveal its other properties. By inquiring into the 
dependence of the mirror symmetry on the angular scale, 
we show that the effect has the highest significance at low 
multipoles i < 5 and that it quickly disappears at higher 
multipoles. 

Similarly, the study of the sum and difference maps 
of temperature inhomogeneity regions provides additional 



insight into the mirroring. Namely, when the sum and 
difference maps are created via reflection of one of the 
maps, as it should be for mirrored images, anisotropic 
properties of excursion sets do survive, while they dis- 
appear if the sum map is created without reflection. 
Difference maps from independent radiometers (A-B) have 
bee n used to test the role of scan inhomogeneities and noise 
(Gu rzadyan et al. 2007b ). The signal-to- noise ratio for the 
studied excursion sets is about 4:1 and the excursion sets 
in (A-B) map do not show any specific property observed 
in the sum (A+B) map. The negligible role of the noise was 
also checked using the foreground reduced maps available 
in http://lambda.gsfc.nasa.gov/product/map/current/, 
Although contamination of a Galactic or interplanetary 
origin at these multipoles certainly cannot be excluded, 
following Gurzadyan et al (2007b), in the last section 
we discuss which properties of the Universe might be 
responsible for this effect, if it had a cosmological origin. 



2. Distribution of excursion sets 

For this analysis we used the 94 GHz (3.2mm) W-band 
WMAP 5-year maps, due to their offering the highest an- 
gular resolution (beam width of FWHM=0°.21) and low- 
est contamination by synchrotron radiation of the Galaxy 
(Bennett et al. 2003j) . The role of the Galactic disk was 
minimized via exclusion of the equatorial belt \b\ < 20°. 

Algorit hms for studying the excursion sets have been de- 
scribed in (Gurzadyan et al. 2 00 5 [[Gurzadyan et al. 2007a ) 
in connection with the study of the ellipticity in excur- 
sion sets in the Boomerang and WMAP, and earlier for 
COBE maps ( ,Gurzadyan fc Torres 1997| ). The definition of 
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the centers and of the geometrical characteristics of the 
excursion sets are based on rigorous procedures, e.g. the 
Cartan theorem for the conjugation of the maximally com- 
pact subgroups of Lie groups. The distribution of the cen- 
ters of the excursion sets obtained via those algorithms were 
obtained for various pixel count and temperature-threshold 
interval sets and the center of these centers was obtained: 
x(center)=^ Xi/N ^ ?/(center)=^ Hi/N . 

Inhomogeneities in the distribution of the excursion sets 
at the temperature interval AT = 90 fiK within |T| = 
45 iiK are concentrated around almost antipodal points 
centered at 

/ = 94°. 7, 6 = 34°.4 {CEn)] 
I = 279°.8, 6 = -29°.2 {CEs). 

At higher temperature interval limits the excursion sets are 
merged and loose their identity. The methodical novelty 
here is that we consider the excursion sets within a tem- 
perature interval [— T, T] with respect to the mean CMB 
temperature, instead of excursion sets in the slice around a 
fixed T. Namely, the symmetry appears at the sum of slices 
but not in each of them. This is analogous to the appear- 
ance of the Great Wall as a result of the sum of slices of 
large scale galaxy surveys. 

Comparing these positions to those of the Maxwellian 
vectors of the lowest multipoles of CMB, it was shown that 
CEjsf and CEs are located close to one of the vectors of 
multipole ^ = 3, shifting towards the equator when increas- 
ing the temperature interval. During this shift, the mirror 
symmetry is approximately maintained, while the patterns 
of the excursion sets around CEn and CEs mirror each 
other with = 0-7 - 1.5. 

Neither CEn and CEs are close to the positions of the 
sum of the multipoles vectors with ^ = 2 — 8, the modulus 
of each vector weighted by l/£{£-\-l)) . The position of CEs 
is not close to that of the cold spot (|Cruz et al. 2QQ5|) . 



3. Mirroring versus multipoles 

We now investigate the multipole dependence of the mir- 
roring. We study it using the method of gradually removing 
multipoles, defined by the coefficients of the tempera- 
ture ffuctuations expansion into spherical harmonics: 



To quantify the degree of mirroring in 
(Gurzadyan et al. 2QQ7b), the drift of the centers of 



aemYemin). 



(1) 



The functions "anafast" and "synfast" of Healpix 
(jGorski et al. 2005^ were used, and the WMAP5 94 GHz 
(W channel) dataset was analyzed. We checked that the 
temperature anisotropy can be represented as 



AT{0, 



AT 



AT 



(2) 



non—mirr 



where the first, mirrored term dominates at certain (low) 
multipoles, while the second, non-mirrored term becomes 
the main one at other (higher) multipoles. Such partial 
mirroring does not imply planarity, i.e. the dominance of 
multipoles with |m| = £. Instead, for the first term on the 
right hand side of Eq. all a^^ with an even value of 
i — m may be non-zero. 



the symmetries with respect to the dipole apexes vs the 
temperatures for each hemisphere was studied. The confi- 
dence levels defined as = (E^^(>"ii - ^2i)^)/(A^ - DOE) 
(where yii,l2i are the distances to the northern and 
southern dipole apexes, respectively, and A/" :^ 30 is the 
number of steps, DOE = 2, for the compared curves) have 
been evaluated, i.e. the departure from being identical (see 
Fig. 4 in ([Gurzadyan et al. 2QQ7bD , as well as Figs. 5 and 
6 there for the distances to the the multipole vectors). 
Here we evaluate the analogeous confidence levels but for 
the maps split over multipoles. Figure 1 shows the plots 
for the excursion sets of more than 50 pixel counts vs the 
multipole numbers where the points of given i indicate 
that multipoles smaller than i are cancelled. The red 
line there represents the best polynomial fit, while the 
dashed lines correspond to smoothed error bars. The plot 
shows that the mirroring effect has its highest significance 



at low multipoles, x < 1.7 for ^ < 5, and it weakens 
monotonically for higher ones; i.e. x^ > 3.5 at i > 10. For 
comparison, we also show the same dependence obtained 
using the WMAP 3-year W-maps (Fig. 2). Similar de- 
pendence occurs for larger excursion sets (Fig. 3). Figure 
4 shows confidence levels for the mirroring effect when 
one of the mirrored regions (the northern, see below) is 
replaced by a simulated, statistically isotropic Gaussian 
map. For the latter, first, the maps for the given multipoles 
were obtained from the real map, then, for the resulted 
multipole maps, isotropic Gaussian maps were generated 
using the mean and the a of each map. 

Then, to probe the mirroring further, we constructed 
and compared the sum and the difference maps of regions 
{l,b: [10°, 170°], [20°, 90°]; [190°, 350°], [-20°, -90°]), re- 
spectively centered on CEn and CEs^ using two proce- 
dures: (a) with rotation on the angle tt, i.e. keeping the 
mirror symmetry; (b) without this rotation. The number 
of the excursion sets surrounding CEn and CEs are given 
in Fig. 5. The angular distances of CEn and CEs from 
the CMB dipole direction vs the temperature threshold are 
shown in Figs. 6a and 5b for the mirrored (i.e. with rota- 
tion) and non- mirrored sums and difference maps, respec- 
tively; the continuous and dashed lines denote the sum and 
difference maps, respectively. In Fig. 6b, the same depen- 
dence is plotted for a simulated statistically isotropic map, 
too (the upper dot-bar curve). Note that, while a difference 
map (A-B) obtained using data from different radiometers 
but from the same region of the sky contains mainly noise 
and no signal, when we deal with different sky regions, the 
difference map should result in another map that is not very 
different from the sum map. Indeed, Figs. 6a and b clearly 
indicate that the similarity in the temperature indepen- 
dence of the distance from the dipole for the non-mirrored 
sum map and the simulated map is obvious - in both maps 
there is no breaking of statistical isotropy. This differs cru- 
cially from the temperature dependence of this distance in 
the case of the mirrored sum and difference maps (Fig. 6a), 
thus confirming the existence of a partial mirror symmetry 
of the regions of CEn and CEs. 

4. Conclusions and discussion 

We have analyzed the scale dependence of partial mirror 
symmetry in the distribution of excursion sets, as previ- 
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Fig. 1. Statistical significance of the mirroring of patterns 
of excursion sets with more than 50 pixel counts around the 
centers CE^ and CEs vs the multipole number i for the 
WMAP 5-year 94 GHz temperature maps. 



ously found in the WMAP3 maps, with a nearly antipodal 
location of symmetry centers. Studies used the WMAP5 
W-band maps, and the results obtained using WMAP5 and 
WMAP3 data agree. The centers lie close to one of the i = 3 
multipole Maxwellian vectors, b ut not close to the sum o f 
multipoles vectors up to £ = 8 ( [Gurzadyan et al. 2QQ7b ). 
Also they are close to the ecliptic pole and are nearly or- 
thogonal to the CMB dipole apexes. They are moving to- 
wards the Galactic equator with the increase in the tem- 
perature threshold interval. 

This symmetry appears to be a large-angle effect, i.e. 
it is stronger at low multipoles and it weakens rapidly for 
larger i: its statistical significance is quantified by < 1-7 
at ^ < 5 and > 3.5 at £ > 10. The symmetry was also 
tested using the following procedure: the sum map of the 
symmetry regions was obtained first, via rotation over tt, as 
is usually the case for mirrored images, and then without 
such a rotation. The clear mirroring in the first case and 
its complete absence in the second case make the case for 
a partial mirror symmetry stronger. 




Fig. 3. dependence as in Fig. 1, but now both vs the 
multipole number and the number of pixel counts of the 
excursion sets. 
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Fig. 4. when one of the symmetry regions is replaced by 
a simulated statistically isotropic map. 
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Fig. 2. The same as in Fig. 1, but for the WMAP 3-year 
maps. 



Turning to the origin of this symmetry, an unknown 
interplay of interplanetary and Galactic foregrounds or an- 
other unspecified non-cosmological contribution to the low 
multipoles certainly remains a possibility. If, however, it has 
a cosmological origin, then a signature of the simplest non 
trivial, T^, spatial topology of the Universe is among the 
options, as discussed in Gurzadyan et al (2007b). For this 
topology, the points with coordinates z and z^L are identi- 
fied where z is one of the spatial coordinates. Such a model 
may be also considered as a limiting case of the Universe 
with compact flat spatial sections having the topology if 
the identification scales Li, L2 along two other spatial coor- 
dinates are much more than L - the slab topology (for early 
papers on a non trivial spatial topology of the Universe, see 
Zeldovich (1973), Sokolov & Schwartsman (1975), Sokolov 
& Starobinsky (1975), Fang & Sato (1983)). Note that, one 
should not expect any mirror symmetry, even a partial one, 
for comparable topological scales ^^^2,^3. 
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As shown in Starobinsky (1993), for this topology, 
a large-angle pattern of a CMB temperature anisotropy 
has just the form The first term on its right hand 
side has the exact mirror symmetry with respect to the 
(x, ?/)-plane. It originates from the Sachs- Wolfe effect at 
the last scattering surface from density perturbations that 
do not depend on z. The second term represents a re- 
maining part of anisotropy and does not have any sym- 
metry at all. However, for a^L on the order of i?hor or 
slightly more, where ao = a(to) is the present scale fac- 
tor of a Friedmann-Robertson- Walker cosmological model, 
the latter term should somehow be suppressed since the 
Sachs- Wolfe contribution to it from the last scattering sur- 
face comes from perturbations having wave vectors with 
|k| > 27r/L. That is why one expects the total large- angle 
pattern of AT/T to have an approximate mirror symmetry 
in this case. E] 

More generically and not connected with a non trivial 
spatial topology of the Universe, such an approximate mir- 
ror symmetry at large angles arises when the large-scale 
^-independent part of density perturbations inside the last 
scattering surface is anomalously large. In all these cases, 
the rms amplitude of the first term on the right hand side 
of Eq. ([2| quickly becomes negligible compared to the sec- 
ond one with the growth of I. Weakening of the mirroring 
symmetry for higher values of I found in Sect. 3 is in a 
good agreement with this theoretical prediction and may 
be considered as an additional argument for the reality of 
the mirroring effect. 

Until now, searches for the mirroring effect of the form 

tor directly for the topology gave negative results 
e e.g. de Oliveira-Costa et al. 1996; de Oliveira-Costa 
et al. 2004; Cornish et al. 2004), rasing the lower limit on 
the physical topological scale a^L to ~ i^hor = 14 Gpc. 
The numerical value is given for the standard ACDM cos- 
mological model with = 0.3, I^a = 0.7. However, for 
higher values of a^L this topology is not excluded. This 



^ Note the other effect worsening the symmetry at large angles 
(low multipoles): a contribution from the integrated Sachs- Wolfe 
effect at small reds hifts due to a cosmolo gical constant (first 
calculated in ( |Kofman &: Starobinsky 1985| )) or dynamical dark 
energy. 
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Fig. 5. The number of excursion sets with more than 100 
pixels vs the temperature threshold within the regions cen- 
tered on CEn and CEs (dashed line) for the WMAP 5-year 
94 GHz temperature maps. 
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Fig. 6. Distances of CEn and CEs from the CMB dipole 
vs temperature for the sum and difference maps centered on 
CEjsj and CEs- (a) The sum and difference (dashed) maps 
obtained via rotation by the angle tt in one of the maps 
(to keep the mirror symmetry); (b) no mirroring rotation 
performed for the sum and difference (dashed) maps; for 
comparison, the case of simulated isotropic maps is shown 
(dot-bar curve). 



follows already from the fact that the much more restric- 
tive cubic topology {L = Li = L2) with a^L > Rhor 
is still considered as a vi able possibility; see the recent 
papers (|Aurich et al. 2007| lAurich 2008|) . where an incon- 
clusive evidence for the latter case with a^L ^ 1.1b Rhor 
is presented. This shows that topological explanation of 
the partial mirror symmetry investigated in this paper is 
possible. From our analysis, it is still too early to speak 
about the value of L, since the secure separation of CMB 
temperature fluctuations into a mirrored and non-mirrored 
parts needs higher resolution maps. Also, a non-topological 
(though still cosmological) explanation of such an effect is 
possible, as pointed out above. In this respe ct, see the re- 
cent papers ( |Gurzadyan fc Kocharyan 2008 ) where it was 
shown that voids can act as hyperbolic lenses in a spatially 
fiat Universe, producing specific signatures in CMB tem- 
perature fluctuations. Future observational data will help 
solve these problems. 

Shortl y before submission of this ma nuscript, a paper 
appeared fCroe neboom fc Eriksen 2008) where CMB sta- 
tistical anisotropy of an axial type was studied with the pre- 
ferred axis very close to the one defined by our CEn — CEs 
direction. 
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